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ABSTRACT 

We derive a viscosity from gravitational instability in self-gravitating accre- 
tion disks, which has the required properties to account for the observed fast 
formation of the first super-massive black holes in highly redshifted quasars and 
for the cosmological evolution of the black hole-mass distribution. 

Subject headings: accretion, accretion disks — hydrodynamics — turbulence - 
galaxies: active — galaxies: quasars: general 

1. Introduction 

Viscous accretion through disks and the ensuing dissipation is known to be a very 
efficient process of converting gravitational energy into radiation. In particular accretion 
into black holes allows to liberate a sizeable fraction of the accreted matter's rest energy. 
For (stationary) accretion at a rate M this amounts to an accretion luminosity L accr of 
-^accr = rjMc 2 . Here c is the speed of light and rj a parameter which takes care of the spin of 
the black hole, i.e., the metrics in the vicinity of the horizon and the corresponding radius of 
the innermost stable orbit (ri SO ), the decoupling of the material from the disk near r iso , and 
the radiation efficiency. For standard accretion disk s, i.e., those which do not belong to the 



class of radiation-inefficient accretion flows (RIAFs: iRees et al. Ill982l ). r] is of order 10 1 
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The evolution of the disk and its timescale are governed by the value of the viscosity 
parameter v. The viscous timescale r visc is given by 

7"vise ■ (1) 

Here s is the radial coordinate in a cylindrical system {R, if, z}. In the following, we will as- 
sume rotational symmetry about the z-axis, and make use of the approximation of vertically 
geometrically thin accretion disks. For exhaustive descriptions of details of the theory of 



this cl ass of disks, we refer the reader to the pertinent literature, e.g-. lFrank. King, fc Raine 



fl20Q2h 



It is not disputed that molecular viscosity is too small by many orders of magnitude 
and leads in almost all relevant situations to timescales surpassing the Hubble time. Far 
less clear is, however , what viscosity to use instead. This impasse was solved originally by 



Shakura fc Sunyaev I (I1973I ). Their ansatz is based on the insight that molecularly viscous 
accretion disks are prone to exceedingly large Reynolds numbers, indicative of the onset of 
turbulence. In their ad hoc prescription, they parameterized the viscous length scale by (a 
fraction) of the disk's thickness h and the velocity by that of the sound speed c s . The entire 
unknown physics was subsumed in a parameter a which was assumed to be (more or less) a 
constant. The assumptions of isotropic sub-sonic turbulence require a < 1. 

This parametrization, v = ahc s , is often referred to as "a viscosity". It proved to be 
very successful for non-self-gravitating disks, i.e., disks in which the gravitational potential is 
solely given by the central accretion body, like in close binaries, late phases of star formation, 
etc. The observed (or derived) evolutionary time scales in these systems allowed to derive 
at least the order of magnitude of a by assuming that the evolutionary time scales can 
be reasonably estimated by r visc . While practically all derived values are compatible with 
the requirement of a being smaller than unity, it also turned out that in the vast majority 
of cases values not too much smaller than this limit were required. Typical values were 
log a = — 1 ± 1 . 

Despite a number of successful applications of the a parametrization (for instance in 
explaining the dwarf nova phenomenon as a disk instability, which is compatible with a 
functional dependence of viscosity on the physical parameters as present in a viscosity), this 
parametrization suffers from a number of shortcomings. 

• It was introduced in a pure ad-hoc fashion, and, in its original form, is not based on 
an instability which could drive the turbulence. While the otherwise exceedingly large 
Reynolds numbers in these flows strongly point towards the occurrence of turbulence, 
the Rayleigh criterion indicates that — at least in the linear regime — the radial angular 
momentum stratification of disks is stable against the onset of turbulence. 
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A naive extrapolation into the regime where the mass of the disk is no longer negli- 
gible leads to the rather unphysical result of a radially constant effective temperature 
(IDuschl. Strittmatter. fc Biermann II2000I ). 



While the energy release from disks is held responsible for the output from galactic cen- 
ters, in particular active ones (AGN), even for an a approaching unity, the timescales 
in these disks are far to long. This problem was stepped up with the detection of 
super-massive black holes (SMBHs) in highly red-shifted quasars. If one assumes that 
these black holes gained their observed masses by accretion, in the most extreme cases, 
less than 10 9 yrs are available for amassing more than 1O 9 M . 



Balbus & Hawlevl (199 



described by IVelikhov I ( 



) re- disco yered a magneto - rotat ional instability (MRI), originally 



19591 ) and IChandrasekhar I (119601 ) . It could serve as an explanation 



for the onset of turbulence, even in only weakly magnetized disks, and led to a viscosity 
of the type and amount proposed by Shakura and Sunyaev. The question whether purely 
hydrodynamic instabilities are also possible in massless disks is not settled so far. While 
Balbus and Hawley's work solved the seeming contradiction between the large Reynolds 
numbers and Rayleigh stability, in particular the importance of non-linear hydrodynamic 
instabilities is far from clear. MRI, however, suffers from the same problems when going into 
the self-gravitating domain. 

The problem with any viscosity prescription of the functional form of a viscosity is that 
there the factor c s is a local quantity, while h, through the vertical hydrostatic equilibrium 
h/R = c s /vtp is actually a quantity that contains global information about the disk's (radial) 
structure: v v is the azimuthal velocity, which in this case is given by Kepler's third law 
(or its relativistic version). Thus, in non- self-gravitating disks, the viscosity prescription 
contains both information about the local and global disk structure. In the self-gravitating 
case, however, the vertical hydrostatic equilibrium depends (almost) only on the local mass 
distribution, i.e., v becomes a purely local quantity. 

The third issue of the far too long viscous time scales in galactic center disks, finally, 
has been tried to overcome by appealing to non-axisymmetric effects, like bars in the disks, 
which are meant to speed up the transfer of angular momentum to larger, and of mass to 
smaller radii by orders of magnitudes. While this, in principle, is a very attractive proposal 
to solve the problem in situations where bars are present, it seems that not in all galactic 
centers bars are present. This then rules it out as a general solution of the problem, though 
bars may very well play an important role in some systems. 



Duschl. Strittmatter. &: Biermann I (119981 . 120001 ) and independently iRichard fc Zahn 



jl999h Droposed a generalization of the Shakura-Sunyaev parametrization which solved the 
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latter two of the above mentioned three major problems of a viscosity, namely the unphys- 
ical transition into the self-gravitating regime, and the exceedingly long viscous time scales 
of galactic center disks. Their ansatz, however, is still a parametrization and not an en- 
compassing solution to th e turbulence problem. Based on laboratory experiments and on 



theoretical considerations (IWendt Ill933l ; iTaylor Ill923l . Il936l ). the viscosity v is written as 



|^ = PRv, (2) 



crit 



where, in analogy with the a parametrization, a scaling quantity (3 = 3?~ it has been defined. 
In this prescription, however, the scaling quantity (3 is not arbitrarily chosen, but is rather the 
inverse of the critical Reynolds number dt CT n, which, in turn, is thought to be of order 10 2 '" 3 . 
In addition to the prescription of eq. (jSJ), it is required that the corresponding turbulent 
velocity scald3 f tur b = vT^V is smaller than or equal to the sound velocity: f tur b < c s . This 
is the so-called dissipation limit. 

The ensuing viscous timescale, r v i SC = with the azimuthal angular frequency 

uo = v^R^ 1 is sufficiently short to allow for efficient disk accretion in galactic centers, thus 
reconciling the observed luminosities and the derived viscous time scales. They are even suffi- 
ciently short to allow for the rapid formation of the SMBHs in the highest red-shift quasars. 
Finally, the transition to non-self-gravitating, dissipation-limited accretion disks (i.e., the 
regime in which Shakura and Sunyaev's original parametrization is so successful) not only 
recovers the a prescription as the limiting case. (3 = 3?~ it ~ 10~ 2 "~ 3 yields corresponding 
values of a which are compatible with the above discussed range of values. 

In this contribution, we discuss gravitational instability as a possible origin for tur- 
bulence, in particular in self-gravitating accretion disks, and its relation to f3 viscosity 
parametrization. In the following Sect., we describe the properties of the instability, and 
in Sect. E] its connection with viscosity. In the final Sect., we summarize and discuss our 
results. 



2. Gravitational Instability in Self-Gravitating Accretion Disks 



For a geometrically thin stationary accretion disk (IFrank. King, fc Raine 1 120021 ) . ne- 
glecting boundary terms, we have 



(3) 



1 The corresponding turbulent length scale is ^turb = \/]3R 
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M 



—2TTRT.VR, 



(4) 



with M the (constant) mass accretion rate, R, tp the radial and azimuthal coordinates, vr/ v 
the corresponding velocity components, and v the kinematic viscosity. Usin g the mono-pole 
approximation for the gravitational potential (jMineshige &: Umemura 1119971 ). approximating 
the disk mass enclosed within a radius R by M&s^R) = R 2 nT l (R), neglecting the central 
mass in comparison to the disk mass, and assuming the /9-parametrization of the viscosity, 
one can express all quantities as explicit functions of M, j3, and R. For the surface density 

2 

(integrated in vertical direction), for instance, one gets S = {^^^j=^J -R -1 - The mass flow 
rate M, in turn, relates directly to the surface density at the disk's outer radius, So = S(-Ro) 
through M = (2vr) 



3 3 
(r/JD2V2 



■^PRqT, 



o- 



I t is well known (e.g. , iMishurov. Peftiev. &: Suc hkov 1976 ;lKato fc Kumar Hl960uKumar 



1960; IStephenson I Il96ll . and in more detail 



Hunter fc Schweiker 



1981 



Hunter fc Horak 



19831 ) that an infinite rotating viscous medium is Jeans unstable, i.e., the Jeans criterion for 
instability 



k < 



2ttGS 



ci 



k 3 



(5) 



is valid, where k is the wave number and G the gravitational constant, kj is the Jeans wave 
number, and Aj := the corresponding Jeans wavelength. 



The same is true for a uniformly rotating viscous disks (ILynden-Bell fc Pringle 



A somewhat more elaborate discussion on viscous disks can be found in 



Gammie 



974). 



19961) 



and lAntonov fc Kondratyev I (119951 ). All this implies that the Toomre criterion for stability 
of an inviscid disk 



Q:- 



C S K 



> 1 



(6) 



(where Q is called the Toomre parameter and k is the epicyclic frequencj|§; iToomre 1964 ; 
Goldreich fc Lynden-Bell I Il965l ) is not a valid criterion for stability in viscous disks. In 



the above mentioned papers, in the Navier-St okes equation the term — — ^2.^5 



AD 



^ is missing, 



which , however, can be of importance (see, e.g.. lSchmit &: Tscharnuter lll995l ; lFridman fc Polyachenko 
1984h . 



Using the Toomre parameter Q as defined in Equation [6], solving for c s and putting this 
into the hydrostatic equilibrium for fully self-gravitating disks, leads to 



Q' 2 



hH 2 



(7) 



2 In the framework of our approximations, n = y/Stt 
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In the framework of our approximations we get: 

Q 2 = 3~. (8) 

This means that, neglecting viscosity, any geometrically thin (h <C R) fully self-gravitating 
accretion disk is Toomre unstable (Q < 1). 

In the following we derive the full dispersion relation for axisymmetric (■§- = 0) waves 
in a thin viscous disk starting from the Navier-Stokes equations in cylindrical coordinates 
with a full handling of the V • a term, where a is the viscous stress tensor (see, e.g., 



Landau fc Lifshitz Ill963h . 



We linearize the disturbances about a stationary value for all quantities X G {£, Vr, v v , $} 
(<3> is the gravitational potential) 

X =X + 5X 

5X = X 1 expi(ut + kR) (9) 
5X <X 

and omit quadratic and higher order terms of small values. Derivatives of small values are 
assumed to be also small. 

In the limit k S> 4, the Poisson equation A5 Q = 47rGS5(z), where 5(z) is the Kronecker 
symbol, is solved by (IBinney fc Tremaine 1119871 ) 

5$ = -2ttG— . (10) 
| A; | 

A linearized combination of the continuity equation for geometrically thin disks (eq. [3]), 
the Poisson equation, and of the axisymmetric version of the Navier-Stokes equation leads 
to the following dispersion relation: 

s 3 + -vk 2 s 2 + ( -v 2 k 4 + 3{l 2 + c 2 s k 2 -2irGi;k)s + isk 2 (c 2 s k 2 + n 2 - 2ttG£A;-3-^U 2 ) = 0, (11) 
3 V3 / R 2 



wher e 9ft(s) = — $$(u). Using the Routh-Hurwitz theorem (see, e.g., iGradshteyn &: Ryzhik 



20001 ) , we determine the sign of the (real) root of [TT] and derive as a necessary criterion for 



stability: 



v 2 

a 3 (k) — elk 2 + Q — 3—zk 2 - 2irGYjk > 0. (12) 

R l 



The only extremum is ko = ^ GY, v i with the two cases 

|^ = 2(c 2 - 3(3 2 vD > =}► c s > 3/3v v or |^ < c s < 3(3v^. (13) 
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For c s > 3f3vtp, ko is a minimum and the condition for stability can be rewritten as 

^ - > L < 14 > 

This condition, however, cannot be fulfilled, since the second term is negative and 3 _1 Q 2 <C 1 
by virtue of eq. [HJ On the other hand, if c s < 3/3v v , then obviously a 3 (k — > oo) < 0, and all 
large k, i.e., small A, are unstable. Thus all geometrically thin FSG /?-disks are unstable. 

A typical dispersion relation is shown in Figure HJ The imaginary part of uo (the second 
solution) is plotted here for the parameters R = 1 ■ 10 18 , Ro = 3 • 10 18 m, So = 30^|, 
c s = 1000— and (3 = 5 ■ 10~ 3 . Solutions 1 and 3 seem to be always stable, i.e., $s(u) > 0. 



3. The Link Between Gravitational Instability and 
Viscosity in Self-Gravitating Accretion Disks 

We propose that the gravitational instability of the disks, instead of a hydrodynamic 
instability is the main driver of turbulence. Neglecting pressure and shear, the wavelength 
where the time scale of viscosity equals the time scale of gravity defines a natural length scale. 
This requires vk 2 = \ / 2nGTik and leads, again in the framework of our approximations, to 

A mi „= (4vr 3 )!/3m (15) 

This characteristic wavelength corresponds, up to a factor of 1.5, to the minimum of the 
imaginary part of u — which is rather independent of all the other parameters — and thus is 
the predominant size of structure in the system. 

If we now identify A m i n with the characteristic length scale of turbulence, we get 

Zturb = v^R = A mi „ = (4vr 3 )^ii?. (16) 
This, in turn, allows us to derive a value for (3 

(5 = JL = 6.5 ■ KT 5 (17) 

Given our approximations, this is a rather remarkable agreement with the value, proposed 
for a hydrodynamically driven turbulence, (3 = 3?~ it = 10~ 2 '"~ 3 . Moreover, one has to note 
that the determination of minimum of the dispersion relation, which of course suffers from 
our approximations, enters by the sixth power. 
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4. Summary 

We have shown (eq. [7]) that geometrically thin, selfgravitating accretion disks are gravi- 
tationally unstable (which, in itself, is of rather little surprise). This instability may lead to 
turbulence and thus viscosity in the disk. This viscosity does not require the presence of a 
(magneto-)hydrodynamic instability. However, although based on different physical consid- 
erations, the functional form of this gravitationally driven viscosity is the same as that of the 
(3 parametrization. For the scaling parameter (3 we derive a value of order 10~ 4 . Thus, this 
viscosity has the required properties to account for the observed fast formation of the first 
super- massive black holes in high ly redshifted quasars and for the cosmological evolution of 



the black hole-mass distribution (IDuschl fc Strittmatter II2005I ). 
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